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Abstract 

The appearances of complex eigenvalues in the spectra of T^T-symmetric quantum- 
mechanical systems are usually associated with a spontaneous breaking of VT. In 
this letter we discuss a family of models for which this phenomenon is also linked 
with an explicit breaking of supersymmetry. Exact level-crossings are located, and 
connections with A/'-fold supersymmetry and quasi-exact solvability in certain spe- 
cial cases are pointed out. 
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1. In a recent paper [0, we discussed the spectra of the followmg family of VT- 
symmetric eigenvalue problems: 



The wavefunction -(/'(x) is required to be square-integrable on the contour C. For M < 2 
this can be the real axis, while for M > 2, it should be distorted down into the complex 
x-plane, as in for the special case a = /(/+!) = 0. 

The principal interest in PT-symmetric problems lies in the fact that, despite not 
being in any obvious sense Hermitian, they often appear to have spectra which are 
entirely real [3-|T2|. In appendix B of a proo/ of this property was given for the class of 
Hamiltonians (||), drawing on ideas from the so-called 'ODE/IM correspondence' p!3|-[l7[| 
(see also [0]). More precisely, it was shown that, for M, a and I real and M > 1, the 
spectrum of (||) is 

• real if a < M + 1 + |2/+1| ; (2) 

• positive if a < M + 1 - |2/+1| . (3) 

Referring to figure |l|, reality was proved for {a, I) G B U C U D, and positivity for 
(a,0 E D. 
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Figure |^: The initial 'phase diagram' at fixed M. 

Special cases of this result were the subject of a number of previous conjectures: for 
a=l=0, M=3/2 in j^; for a=l=0 in Q; and for a=0, I small in |]l^. More unexpected 
was the possible loss of reality in the region A. The condition (|2|) arose for technical 
reasons in the proof in though it was verified that at some (but not all) points inside 
A, reality did indeed break down. But the reason for this breakdown, and the physical 
significance - if any - of the boundaries of the region A, were left obscure. 

In this note we show that the appearances of complex eigenvalues as the lines 
a = M-|-1-|-|2Z-|-1| are crossed are directly related to the fact that on these lines 
the problem (||) can be reformulated in terms of a T^T-symmetric version of supersym- 
metric quantum mechanics, and exhibits level-crossing. Elsewhere in the phase diagram, 
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supersymmetry is explicitly broken and the level-crossing is lifted. Moving into the re- 
gion A, this lifting occurs through the creation of a pair of complex-conjugate energy 
levels. Since the presence of complex energy levels indicates the spontaneous breaking 
of VT symmetry we see that the model shows an interesting interplay between the 
two different ways in which a symmetry can be broken. 

2. To streamline some formulae, we replace x by x/i in (|^) by setting ^{x) = ij{x/i), 
so that the eigenproblem becomes 



dx"^ 



= -A«>(x), ^{x)eL\iC). 



(4) 



We now specify the quantization contour more precisely: iC starts and ends at |x| = oo, 
joining the (Stokes) sectors 5_i and Si, where 



arg(x) 



27r/c 
2M+2 



This is illustrated in figure 
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Figure]^: The quantization contour iC (arrowed line). 

It will also be convenient to adopt a new set of coordinates on the (a, I) plane, by setting 

1 



2M+2 



[a-M -1± (21+1)' 



(6) 
(7) 



so that 

a = {M+l){l + a+ + a^) , 21 + 1 = {M+l){a+ - a.) . 

The domain boundaries - the dotted lines on figure || - are the lines a± = 0. On the 
line a+ = 0, a = M—21 and the problem (^) factorises as 



Q+Q_$(x) = -A$(x) , with Q± 



±— + x 
L dx 



M 



(8) 



Such a factorisation is usually taken to signal a relationship with supersymmetry |p^ , |20| . 
Indeed, it is immediately seen that (0) has a A = eigenfunction in L^(iC) : 



■^{x) =x 'exp {jjj^: 



,M+l 



Q_-^{x) = 0, 



(9) 
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which can be interpreted as having unbroken supersymmetry. All other eigenfunctions 
are paired with those of the SUSY partner Hamiltonian 7i = Q_ Q_|_ . This is found 
by replacing {a, I) = {M-2l,l) by (aj) = {-M-2l,l-l), or (a+,a_) = 
by (q!+,3_) = (— 1, ~1~ = (— 1, Q- — (A rather different point of view 

on supersymmetry in T^T-symmetric quantum mechanics was taken in where the 
action of the supersymmetry generators Q± was supplemented by T, resulting in a set- 
up with two states of zero energy, rather than one.) Replacing — / by /+1 gives a A = 
eigenfunction on the line a_ = 0. Thus the boundaries between the regions on figure |^ 
are picked out by the presence of a (supersymmetric) zero-energy state in the spectrum 
of the model. One would normally expect this to be the ground state, and indeed this 
is the case on the boundary of D. However, level-crossing means that ^{x) is only 
the ground state on the boundary of A for a < M+3. There is no contradiction with 
the usual theorems of supersymmetric quantum mechanics, since the problem under 
discussion is not Hermitian. 

To verify that level-crossing does occur, we consider T(0), where T(— A) is a spectral 
determinant which vanishes if and only if (^) has a solution, square-integrable on iC, at 
that value of A. In an expression for r(0) was found. In the 'light-cone' coordinates 
a± , this is 

As expected from supersymmetry, T\x=o is identically zero when either a_|_ or a_ van- 
ishes, on account of the presence of the state ^ in the spectrum. Level-crossing will 
occur when a further level passes through zero, but the presence of ^ makes this hard 
to detect from an examination of T(— A, a+, a-)\\=Q alone. It is tempting suppose that 
the level-crossings happen at the double zeroes of T|a=o- Tempting, but wrong, for 
reasons which may become clearer when figures ^ and |^ below are examined. Better 
is to consider the SUSY partner potential to (^), which is isospectral to it save for the 
elimination of the state ^. Substituting the values of S+ and a_ into (|lO|), 



M+1 



a- 



Level-crossings with the state ^ are indicated by simple zeroes of (11), and are at 



(a+,a_) = (0,n + f^) , n = 0,l,.... (12) 



(Had we looked instead for double zeroes of (|lO| ) , we would have predicted - incorrectly — 
the values a_ = n .) Swapping Q!-|_ and a_ throughout gives the level-crossings on the 
line a_ = 0. Note that these level-crossings are exact - the state ^' is protected by 
supersymmetry, and cannot mix with any other state, even as the level-crossing value 
of / is approached. However, as soon as the supersymmetric line is left, the protection 
is lost and mixing does occur. This is the key point, showing why the boundaries of the 
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region A are of more than just technical significance. Since the model is "PT-symmetric, 
all energies are either real or occur in complex-conjugate pairs Q, and a real energy can 
only become complex if it first pairs off with another real energy. What we have just 
shown is that the hidden supersymmetry of the theory on the lines a+ = and a_ = 
affords a mechanism for this pairing-off to occur, by permitting eigenvalues to become 
exactly-degenerate. 

3. To verify this picture, we report some numerical data for M = 3. This value is chosen 
principally for convenience, since at M = 3 the fifth spectral equivalence discussed in ^ 
allows the 'lateral' eigenvalue problem (Q), defined on the contour iC, to be mapped 
onto a radial problem defined on the positive real axis. (Note that this mapping does 
not rule out the appearance of complex eigenvalues of (^), since at such points the 
corresponding radial problem has 'irregular' boundary conditions at the origin and is 
not self-adjoint.) The merit of the radial problem is that it is straightforward to treat 
numerically, using (for example) the Maple code in appendix A of j^. 
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Figure |^: The appearance of complex eigenvalues as the supersymmetric 
line a+ = is crossed, for M = 3. The eigenvalues are plotted against a_ , 
with the points a_ G Z+, E = indicated by boxes; for — 0, these are 
double zeroes of r(— A, a+, a-) . 

In figure ^ we show the real energy levels along the supersymmetric line q_(_ = 0, 
and along two lines just either side of it. On the third of these plots - the only one 
partially inside the region A - pairs of eigenvalues join and become complex in the 
neighbourhoods of the level-crossing points. By contrast, on the first plot, which lies 
entirely inside regions covered by the proof of the spectrum remains entirely real 
and the effect of super symmetry-breaking is seen in the replacement of level-crossing by 
level-repulsion. 

It might seem surprising that we are observing level-crossing in such a simple 
quantum-mechanical system. Again, the non-Hermiticity of the problem provides the 
explanation. As the level-crossing point is approached, not only the two eigenvalues 
but also their eigenfunctions become equal. (This is clear since eigenfunctions of (^) are 
uniquely characterised by the values of a± and E, and their decay in the sector 5_i.) 
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The Hamiltonian thus ceases to be completely diagonalisable, a situation which can be 
understood by considering the following 2x2 matrix, with eigenvalues ±V6 : 



1+6 
6 



1 



±V6 



(13) 



The value of 6 corresponds to the distance from the level-crossing in the supersymmetry- 
breaking direction. (The addition of diag(— r/, 0) to Ti-^^Q would instead model a per- 
turbation which preserves supersymmetry, leaving the theory on the line a+ = 0.) At 
6 = 7] = the two eigenvectors coincide. For this reason, the reference at the end of the 
last section to a degeneracy in the model at the level-crossing should be treated with 
caution - better maybe to say that the eigenproblem itself has become singular. The 
two eigenvalues become complex for 6 negative, so the transition to complex eigenvalues 
can be traced to the singular nature of the eigenproblem at the level-crossing. 

(To avoid confusion, we should stress that we are talking only of the 'bosonic' part 
of a fully supersymmetric problem here - the full SUSY QM system involves both 
Ti.{a+, a_) and "^(3+, a_), and the singularity of the eigenproblem at the level-crossing 
is reflected in the fact that the A = degeneracy jumps from 1 to 2 there, rather than 
from 1 to 3 as would naively have been expected.) 

The lines a± = are not the only ones along which there is a 'protected' zero-energy 



level in the spectrum of the system - from the formula (10), the same is true of all the 
lines a± = n, n £ Z^. At least for M = 3, we can understand this as being due to a 
hidden AA-fold supersymmetry in the model [||,|2^. In figure ^ we illustrate how further 
complex levels are created as these lines are crossed. 



a+ = 0.99 a+ = l a+ = 1.01 

Figure ^ Yet more complex eigenvalues for M = 3, this time near the 
line q;+ = 1 . Labelling as in figure |^. 

Three of the low-lying levels in the middle plot of figure two of which become complex 
for a_ > 1/4, can in fact be found exactly: they are the roots of 

P3{E) = E^ - 32[(5 - 3J)a_ + J - 2]E 







(14) 



at J = 3. More generally for M=3, at a+ = (J— 1)/2, J = 1,2, ... , J levels of the 
"PT-symmetric problem (^) can be found exactly. This can be seen using the above- 
mentioned fifth spectral equivalence of [ffl]: on the lines a+ = (J— 1)/2, (13) is mapped 
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to the radial problem at a point where it is quasi-exactly solvable [p2U[*|. (Equivalently, 
it follows directly from the termination of the relevant Bender-Dunne expansion.) 
For J odd, one of these QES levels is the exactly-zero level mentioned in the previous 
paragraph. This level (together with the other J— 1 QES levels) can be eliminated by 
passing from (a+,a_) = ((J— l)/2,a_) to (a+,a_) = (— (J+l)/2, a_— J/2) [||, and 
the exact locations of the level-crossings can again be found: they are at Q!_ = J/2 + n, 
n = 0, 1, ... . However, for the moment we do not know if a similar exact treatment 
can be given for other values of M, though we expect the qualitative features of the 
breaking of level-crossings to persist, at least for nearby values of M. (Recall that the 
'protection' of the zero-energy state for a± G Z+ holds for all M.) 

To end this section, we show three further plots which should help the reader to 
understand how the levels reorganise passing between figures ^ and ^ The middle plot 
is another QES example. This middle plot illustrates an additional feature of the models 
on the QES lines M=3, a± = (J— 1)/2 : the levels which become complex along these 
lines always lie in the QES part of the spectrum. To prove this, we simply note that 
the 'dual' problems, with {a±,a^) = {a±,azip) = (— (J+l)/2, a_— J/2), always lie in 
regions of the parameter space covered by the reality proof of As already mentioned, 
the spectrum of this dual problem is identical to that of the original problem, minus 
the QES levels. (Strictly speaking, the arguments of ^ applied to a radial problem, 
but the discussion of section 8 of that paper, based on intertwining operators, can be 
modified to cover the current, lateral, problem. One can alternatively argue via the 
'fifth spectral equivalence' mentioned at the start of this section.) This demonstrates 
that the non-QES levels always remain real. A similar result was conjectured for quartic 



QES potentials in |25|, but so far as we know remains unproven. 
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Figure |[ The reorganisation of the levels passing from a+ = to 
a+ = 1 . Labelling as in figure ||. 



4. The above discussion has confirmed that the domains where the spectrum of (Q) 
has a complex component open out from the level-crossing points on the boundary of 

* We remark that a relationship between generalised supersymmetry and quasi-exact solvability, albeit 
for a different set of models, was originally discussed in |p3]. 
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the region A. These points are given exactly, for general M, by equation (12). It 
is interesting to see, at a qualitative level, how these 'domains of unreality' coalesce 
as one moves further into A. Figure ^ below exhibits some initial numerical results, 
found again for M = 3 using the Maple code of Q. This should be viewed as a 
refinement of the initial 'phase diagram' of figure |l[ The full domain of unreality is the 
interior of the curved line, a proper subset of A which only touches its boundary at 
the points {a±,a^) = (0,n+l/2), which are the level-crossing points (^) for M=3. In 
the small, approximately-triangular region inside A but outside the curved line which 
abuts the point q+ = a_ = 0, the spectrum is not only real but also entirely positive, 
despite the fact that it lies outside the domain D. This shows that, while the condition 
a < M -\- 1 — |2/-|-1| is sufficient for positivity of the spectrum, it is not necessary. 

Perhaps the most striking feature of figure ^ is the pattern of cusps. In the absence 
of an analytical analysis, we do not yet know whether these are special to M = 3, or 
more generic. At least numerically, they lie exactly on the lines a± = n, along which the 
model possesses a protected zero-energy state (to guide the eye, segments of these lines 
have also been added on figure |6|). Furthermore, for M = 3 the model is quasi-exactly 
solvable on these lines. As shown above, the levels which go complex then lie in the 
QES part of the spectrum, and this allows the cusps to be located exactly for M=3, 
assuming that they do indeed lie on the QES lines. For example, this places the two 
lowest pairs of cusps at (a, I) = (9, -l/2±3/2) and (15-3/\/2, -l/2±3/\/2) . 




Figure |6[ The domain of unreality for M = 3. Horizontal 
and vertical axes are / and a respectively, as in figure |l|. 

A detailed study of these properties must await future work; we expect that complex 
WKB techniques will be useful in this regard. These were recently employed in |^] for 
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a problem which can be mapped onto the I = 0, M = 2 case of (Q), save for a different 
choice of Stokes sectors - 5_i and ^2 - for the quantisation contour. 

The intersections of the hues 0+ = (Ji — 1)/2 and a_ = (J2— 1)/2 provoke one further 
thought, for M=3. At these points the model is quasi-exactly solvable from two different 
points of view, with either Ji or J2 levels lying in the QES part of the spectrum. This 



is reflected in a curious factorisation property of the Bender-Dunne [24| polynomials 
Pj{(^±iE). (These polynomials encode the QES levels in their zeroes, equation (14) 
being one example.) Taking Ji < J2, at the point (a_|_,a_) = ((Ji— 1)/2, (J2— 1)/2) 
where the two QES lines intersect, Pj-^{a-, E) is a factor of Pj^ia^, E) . A direct proof 
of this result, using the Bender-Dunne three-term recursion relation |2^, can also be 
given. 

5. We conclude with some general comments. Our main purpose in this note has 
been to understand the physical reasons for the breakdown in the reality proof of [Q]. 
We have shown that, on the lines along which the proof fails, the model has a hidden 
supersymmetry, and that once these lines are crossed, the reality property can, and at 
some points does, fail to hold. This implies that the conditions required by the proof 
are of more than technical significance. It would be interesting to extend the proof to 
cover the full domain for which the spectrum of (jl]) is real, but the complicated shape 
of the domain of unreality shown in figure ^ suggests that this will be a difficult task. 

The model we have been discussing has turned out to have an unexpectedly rich 
structure, and should serve as a testing-ground for other aspects of PT-symmetric 
quantum mechanics. More detailed studies should help us to gain a better understanding 
of the emergence of complex eigenvalues in general, as well as shedding light on the 
pattern of transitions revealed by figure |6[ We should also remember that, via the 
ODE /IM correspondence []l3|-[l7[| , all of these results have potential implications in the 
field of integrable models in 1+1 dimensions. 

Finally, we remark that similar appearances of complex levels preceded by level- 
crossings have been observed in a nonunitary model of quantum field theory - the 



boundary scaling Lee- Yang model [26|. In this case the level-crossings were protected 
by an identity between cylinder partition functions in models with differing boundary 
conditions, which followed from a set of functional equations called a T-system. It 
remains to be seen whether this phenomenon can also be understood on the basis of 
some hidden symmetry, as was the case for the level-crossings discussed in a simpler, 
quantum-mechanical context in this letter. 
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